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Abstract 

We study higher-dimensional neutrino mass operators in a low energy theory that contains 
a second Higgs doublet, the two Higgs doublet model. The operators are relevant to underlying 
theories in which the lowest dimension-five mass operators would not be induced. We list the 
independent operators with dimension up to nine with the help of Young tableau. Also listed are 
the lowest dimension-seven operators that involve gauge bosons and violate the lepton number 
by two units. We briefly mention some of possible phenomenological implications. 
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1 Introduction 



The tiny neutrino mass can be accommodated at low energies by nonrenormalizable, higher- 
dimensional mass operators. With the lepton fields as established in the standard model (SM) 
and the Higgs fields assumed to be a doublet, such operators first appear at dimension five flU: 

^=^H*aH f p F Ly , &f=^lfilH a . (1) 

Here H a is the a-th Higgs doublet with hypercharge Y = +1, and Fi x is the %-th left-handed 
leptonic doublet with Y = — 1. A tilde denotes the complex-conjugated field that transforms 
under SU(2)i exactly as the original one, e.g., Fi = io 2 F[, while the superscript C denotes 
charge conjugation with the convention F^ = (Fl) c . 

Both operators G and break the lepton number by two units. When the neutral compo- 
nents of the scalar doublets develop a vacuum expectation value (VEV), & generates a mass 
for neutrinos that is inversely proportional to the energy scale A of some underlying theory 
responsible for the operator. Although the operator & does not generate a mass but involves 
interactions amongst leptons and scalars of different charge, it may arise from the same mech- 
anism that induces 6 due to the similar structure. With a single Higgs doublet as in SM, the 
operator & is unique while 0? does not exist. 

It is interesting to realize that the unique operator 6 in SM may be written in three appar- 
ently different ways 0. This amounts to forming a singlet in three ways out of four factors of 
the two half-isospin fields, and suggests its possible origin from three types of seesaw mech- 
anisms 01 SI [5J. A phenomenological issue with those mechanisms is that the energy scale A 
is so high that it would not be possible to detect any other effects pertained to the origin of 
neutrino mass. From the viewpoint of effective field theory, the scale may be lowered if the 
mass is induced not from a dimension-five operator but from those of even higher dimensions. 
It is conceivable that there will be more and more mechanisms that can induce a mass oper- 
ator as its dimension increases, see Refs. |[6l 13 for some recent examples. However, it has 
been established recently that the mass operator at each higher dimension is always unique [H. 
This implies that as far as the neutrino mass is concerned different mechanisms are completely 
equivalent. But with a lowered scale, it becomes possible to distinguish them through other 
effects. 

In this work we will address the neutrino mass operators in an effective field theory that 
contains two Higgs doublets. Although the two Higgs doublet model (2HDM) is interesting in 
itself, the main motivation comes from supersymmetry which is a leading candidate for physics 
beyond SM and is under examination at high energy colliders. It would be also tempting to see 
how those higher-dimensional operators are induced in a supersymmetric framework. We will 
show that with two Higgs doublets the operators are no more unique but increase quickly in 
number with their dimension. We will list all mass operators of dimension up to nine as well as 
related dimension-seven operators involving the SM gauge fields. We will also discuss briefly 
some of the phenomenological implications of these operators at low energies. 
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2 Mass operators up to dimension nine 



We assume that the low energy theory contains the SM fields and an additional Higgs doublet 
that also develops a VEV. The neutrinos can only have a Majorana mass in this case. We are 
interested in the high-dimensional operators that can yield a neutrino mass (called i^-type) when 
the Higgs fields assume their VEV's, as well as those that do not give a mass but have a similar 
structure (^-type). We can therefore restrict ourselves to the two-lepton sector that violates 
the lepton number by two units. The relevant fields are the lepton doublet Fl, the two Higgs 
doublets H\ and H 2 plus their properly complex-conjugated fields which also transform as a 
doublet under SU (2)/,: 



a = F£(-l),b = F L (-l), c = H x (+1), d = H 2 (-l),e = H 2 (+l), f = H l (-1), (2) 



where the number in parentheses indicates hypercharge. Our notation is such that we always 
use column spinors in isospin space though F^ is a row spinor in Dirac space and should appear 
on the left of Fi to form an appropriate Dirac bilinear. The lepton generation index is generally 
inessential and can be easily recovered when necessary. We note the following features that 
are useful to exhaust all possibilities. First, since the pair ab appears once, there are two more 
factors of c or e than d or / to balance hypercharge. The dimension of mass operators is thus 
In + 5, where n denotes the number of copies of d or /. Second, the occurrence of c may 
be replaced by e if this yields a different and nonvanishing result, and similarly with d and /. 
Finally, the SM case is recovered by the identifications e = c and d = f. 

With an even number of fields with nonzero isospin one may imagine to form higher isospin 
products before building a singlet out of them. But this is unnecessary when all the fields are in 
the fundamental representation (spinor for short) of SU (2) : all isospin invariants of a given mass 
dimension can be exhausted by first forming singlets from any two spinors and then multiplying 
them. This is the group-theoretical reason that the three types of seesaws reduce to the unique 
dimension-five Weinberg operator 6 in SM and that its higher-dimensional generalizations 
are also unique at each dimension [8l . 

The above point can be best seen in the tensor method in terms of Young tableau. For 
SU (2) a Young tableau has at most two rows, and each column with two rows is a separate 
invariant. This is especially convenient when only spinors appear, because in that case each box 
represents an individual field and a two-row column is an antisymmetric, invariant product of 
the two spinors involved. This has a few immediate consequences. First, there can be no bare 
mass term from F^F^ even if F^ had a zero hypercharge. Second, denoting a spinor by its index 
in the box, we have the basic relation: 



i m 



LB = ±1 



n\j\ \m\n 



(3) 



which is equivalent to the relation (z , j, m, n = 1 , 2) 




(4) 



(5) 
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which means that only one group of dimension-five operators (type 0) listed in Ref [QQ| are 
actually independent. (Be careful not to mix the generation indices with the spinor indices.) 
More generally, putting spinors directly in boxes we have 



a 


K 


b 


a 



a 


b 


K 


<7 



a 


b 


a 


K 



(6) 

where W is any Young tableau. Namely, the ^-type operators that contain as a factor an 
invariant formed out of a, b are linear compositions of the ^-type operators. By making a 
complete list of all mass operators (of type-^), all non-mass operators (of type-<5^) with a 
similar structure are automatically covered. In the language of Young tableau, we will never put 
a, b in the same column. 

It is easy to figure out all dimension-five operators since J, / cannot appear while c/e 
appears twice. They are 



a 


b 


c 


c 



a 


b 


c 


e 



(7) 



plus those obtained by c -H- e, or 



S 5 = (a,c) (b,c)o, T 5 = (a,c) (b,e) , 

$5 = ^Icoe? 2fi = ^5|e«-fej 



(8) 



where the subscript denotes an isospin invariant formed by antisymmetrizing the fields inside 
the parentheses; for instance, denoting the upper (lower) component of a spinor by a subscript 
plus (minus) sign, we have y/l(a,c)o = a + C- —a^c + . Since a and b are essentially the same 
field, the list of operators may be further reduced. To see this clearly, we reserve the lepton 
generation index by putting a = F^ x , b = Fi y . Then, 



= ~ f^+) (vLyC- - f Ly c + ) = {vgc- - fgc + ) (v^_ - /z^+) = 22? 



yx 



(9) 



where V^Vj = ls used. We can thus choose 5s, §5, T$ as the complete and independent 
list of dimension-five operators. 

At dimension seven, the operators contain three copies of c or e and one copy of d or /, 
and can be classified as S : c 3 d, T : c 2 ed, plus those obtained by c <H> e, or d <H> /, or both 
interchanges. The first one is easy to write down: 



S 7 = (a,c)o(b,c)o(d,c)o. 



(10) 



For the second one, there are following possibilities to distribute the spinors in the boxes of a 
2x3 Young tableau: 



a 


b 


J 


c 


c 


e 
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b 
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c 


e 


c 



st 



2nd 



a 


b 
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c 


d 


c 



3rd 



a 


b 


d 


e 


c 


c 



4th 



a 


b 


e 


d 


c 


c 



5th 



But the basic relation in eq © implies 

lst-2nd + 3rd = 0, 1st- 4th + 5th = 0, 



(ID 



(12) 



4 



which eliminate two operators. We choose the 1st, 3rd, and 5th ones to be independent: 

T 7 l = (a,c)o(b,c) (d,e)o, T 2 = (a,c)o(b,d) (e,c) , T 3 = (a,d) (b,c)o(e,c) . (13) 

But for the same reason as for T$, T 3 is covered by T 2 when the lepton generation indices are 
reserved, and may thus be excluded as redundant. The remaining operators are obtained by 
interchanges: 

£7 = SVlc-H-ej £7 = S^U-H-/) ^7 = SVlcOe^-B-/; 

'f'l^ yl,2i 7^1,2 t^'^I '7^t)2 t ' \ C\A\ 

There are altogether 12 operators at dimension seven. 

The dimension-nine operators contain four copies of c or e and two copies of d or /, which 
are classified as S : c A d 2 , T : c 4 df, U : c 3 ed 2 , V : c 3 edf, W : c 2 e 2 d 2 , X : c 2 e 2 df, plus those 
obtained by interchange c -H- e, or J -h- /, or both. We continue to denote an operator obtained 
by c <H> e with a bar, that by d <H> / with a hat, and the one by both c <H> e and J <H> / with a tilde. 
It is easy to write down S and T: 

S 9 = S 7 (d,c) , T 9 =S 7 (d,c) . (15) 

And there are four more operators that are obtained by interchanges: 

S9, S9, 5g; Tg. (16) 

The U operators have one more factor of cd than Tj. Deleting the redundant one associated with 
T 3 , we have 

U g l ' 2 = T 7 l ' 2 (d,c) . (17) 
The additional operators obtained by interchanges are also independent 

U g l > 2 , ft}' 2 , U 1 /. (18) 

It is possible to distribute in nine ways the spinors of V in a 2 x 4 Young tableau, but only 
four of them yield independent operators. Five Young tableaux are obtained from those in eq 
(ITT]) by attaching an additional column of (/, c)o, thus giving the three independent operators: 

^1,2,3 = z ,l ) 2,3 (/jc)o _ (19) 

Another three tableaux are obtained from the above by d -H- /: Vg 1 ' ' , and the ninth one corre- 
sponds to 

Vg = (a,c)o(b,c)o(d,f)o(e,c) . (20) 

Similar to T? and T 3 , Vg can be deleted from the list. Furthermore, the basic relation © implies 
that 

Vg 1 - Vg 1 = Vg 2 ~ Vg 2 = Vg 3 - Vg 3 = -Vg , (21) 
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12 ^ 1 2 3 

so that we can keep Vg ' ' in the list while excluding V 9 ' ' as redundant. Finally, there are three 
more operators obtained from interchange c -H- e: 

Vg' 1 ' 2 . (22) 

A similar (but slightly different) analysis applies to the operators W, which have four indepen- 
dent forms 

W 9 l ' 2 = T 7 1 ' 2 (d,e) , W 9 1 ' 2 , (23) 
plus four more by interchange d <H> /: 

Wg' 2 , W 9 > 2 . (24) 

Finally we come to the symmetric case of X that has the most possible Young tableaux (18 
in total). The basic relation © removes ten of them as redundant and the symmetry in the 
lepton fields deletes another three, leaving us with five independent operators: 

Xg 2 = T^ 2 (f,e)o, Xg 

Xg = (a,c)o(b,e) (d,f)o(c,e)o, 

Xg 5 = (a, d)o(b, /)o(c, e) (c, e) , (25) 

where Xg (Xg) is (anti)symmetric in c -H- e and d f respectively when the lepton generation 
indices are ignored. There are altogether 33 dimension-nine mass operators. 



3 Adding gauge bosons 

The underlying physics that produces the higher-dimensional neutrino mass operators in the last 
section may also induce lepton-number violating interactions with gauge bosons. In this section 
we continue to work in the two-lepton sector and list the lowest dimension-seven operators with 
gauge bosons that are built upon the dimension-five mass operators. The gauge fields may 
enter in two ways, either through gauge covariant derivatives or through field strength tensors. 
The first case amounts to introducing new Lorentz vector fields that have the same quantum 
numbers under the SM gauge group as the original fields, a, b, c, e. The second case requires 
that those original fields must be built into a hypercharge-neutral, i so spin-triplet or -singlet form 
that couples to the field strength tensors of SU (2)l and U(1)y respectively. 
We start with the operators containing the gauge covariant derivative 

£> M = d^-ig 2 ^a a W^i 8l ^B^ (26) 

where the minus (plus) sign applies to the fields c, e (a, b), and Wfi and are the gauge fields 
with gauge couplings g2,l- Distributing two factors of to any two of the four fields in the 
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mass operators S5 and T$ yields 



7 1 '-' 6 = (D^a,D^c)o(b,c) , (D^a,c) (D^b,c) , (D^a,c) (b,D^c) , 

(a,D^c) (b,D^c) , (a,D^c) (D^b,c) , (a,c) (D^b,D^c) ; (27) 
K 1 '-' 6 = (D^a,D^c) (b,e) , (D^a,c) (D^b,e) , (D^a,c) (b,D^e) , 

(a,Z) M c)o(&,Z)^e)o, (a,D M c) (^,e)o, («, c) (D M &,D M e) ; (28) 

plus J 1 ' "' 6 and £ 1 ' "*' that are obtained by c -H- e. Since the gauge covariant derivative does 
not spoil the relation ty/f ty/; = tyftyi, we can exclude some of the operators as redundant as we 
did with 75 . Reserving the lepton generation indices and denoting the upper (lower) component 
of a gauge covariant derivative also by a subscript plus (minus) sign, we have, for instance, 
a x+ (D^b y )- = (D^a y )-b x + using our notations in eq ©. It should be reminded that no inte- 
gration by parts can be legitimately used here; instead, the relation y/f 1//)" = V^f tyi 1S sufficient. 
Some inspection then shows that = J y l x , J~y = Jy X and similarly for /. Since the ^operators in- 
volve simultaneously c and e fields, a stronger reduction of operators becomes possible, namely, 
Kxy ' ' ' ' = K yx ' ' ' ' . The complete and independent operators can thus be chosen to be 

7 1 '"' 4 , J 1 ''"' 4 , K 1 '"'* 6 . (29) 

To construct dimension- seven operators involving gauge field strength tensors, the Lorentz 
indices of the tensors must be contracted by Dirac matrices. This means that a o^ v should be 
sandwiched between the lepton fields a and b, which fits well with their chiralities. Consider 
first the coupling of B^ v to a singlet formed from abce. The only difference to the dimension- 
five mass operator T$ is to insert a a^ v between a and b: 

M(B) = (a,c) o> iV (b,e)oB ll y. (30) 

Note that M(B), which is again obtained from M(B) by the interchange c <H> e, is not indepen- 
dent since when attaching the lepton generation indices to a and b we have M^B) = —M yx (B) 
upon using i/^cfyvV'y = ~¥y°^vVx- It is not necessary either to consider the case where a, b 
lie in the same column of a tableau since the basic relation © is not disturbed by the Lorentz 
structure. Similarly, the counterparts of S$ and S5 are 

L(B) = (a,c)oa^ v (^,c)o J S Alv , L{B). (31) 

In contrast to the above, the SU (2) L gauge field strength W£ v being a triplet must couple to 
the triplet states of abce to become a singlet. There are apparently nine possibilities for abce to 
form a triplet state: one pair of spinors in a singlet and the other in a triplet (six in total), or both 
pairs in a triplet multiplied into a triplet (three in total). But only three of them are independent 
as we show below. We note first of all that there are four possible ways to form a state with 
/3 = +1, which do not necessarily have a definite total isospin /: 

w = a_b + c + e + , x = a + b_c + e + , y = a + b + C-e +1 z = a + b + c + e_. (32) 
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But symmetry requires that the 1$ = +1 state of a triplet (1=1) formed from four spinors be 
a difference of the above quantities, and thus there can only be three independent states with 
/3 = +1 belonging to three triplets. (The fourth I^ = +1 state has 1 = 2, and all this is consistent 
with isospin composition indeed.) For instance, the Ij, — +1 state formed from ab in a singlet 
and ce in a triplet is (w — x)/ y/2, while the I3 = +1 state with all of ab, ce, and abce in a triplet 
is given by (w + x — y — z)/2. 

To write isospin- 1 states formed with two isospin-half ones, it is convenient to use the row 
spinor. We denote by a check the combined action of tilde and dagger on the isospin space, 
which transforms a column spinor to a row spinor in the complex representation. For instance, 

b = P = (f L ,-v L ), (33) 

can form a singlet with c and the gauge field strength, bW^ v c, where W^ v = j(J a W£ v . The 
complete and independent couplings of W^ v to abce are therefore as follows: 

M l (W) = (a,a^ v b) cW^ v e, M 2 (W) = {a,c) Q o llv bW ilv e, M 2 (W). (34) 

M 2 (W) is independent of M 2 (W) since c and e are now at inequivalent places in contrast to 
M(B). On the other hand, M l (W) = M l {W) because cW^e = eW^c. This is in accord with 
the above symmetry arguments. The operators with two c or two e are 

L\W) = (a,o^ v b) cW^ v c, L 2 (W) = (a,c) o^ v bW^ v c, L l (W), L 2 (W). (35) 

To summarize, the complete and independent dimension-seven operators involving the gauge 
field tensors are L(B), L(B), M(B), L l (W), L 2 (W), L l (W), L 2 (W), M 1 (W), M 2 (W) and 
M 2 (W). 



4 Discussion 

The effective operators that we have written down in the last two sections involve various 
lepton-number violating interactions of multi-Higgs and gauge bosons, which may have rich 
phenomenological implications. But to make a complete analysis, we should include some 
other operators at a similar dimension, in particular those involving four-fermions, that violate 
the lepton number by two units. Operators involving four and six fermions in SM were ana- 
lyzed in Refs. |[9l [TOj for inducing neutrino mass at the loop level and their phenomenology 
explored in ifTOl . The neutrino mass operators with two Higgs doublets were symbolically writ- 
ten down in Ref. from hypercharge balance but no attempt was made to complete their 
isospin structures. Instead, possible underlying models were suggested that could induce a 
specific dimension-seven operator via seesaw, together with radiative mechanisms. 

In this concluding section, we discuss briefly some interesting interactions contained in the 
operators listed in the last sections, while leaving a more complete phenomenological analysis 
for the future work. Assume both Higgs doublets develop VEV's which are generally complex 
with phases «i ; 2, 

(flf > = ^=u lCp , (£#) = ^=u 2 s p , (36) 
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where v = 246 GeV and cp = cos/3, sp = sin/3. The would-be Goldstone bosons G^ and 
physical scalars H ± 1 A , R a (a = 1, 2) are related to the original fields by unitary transforma- 
tions: 



H ) \ u\cp u 2 sp J \ H{ 
il a \ 1 / —Ua \ ( 



c j8 S P J \ h 



(37) 
(38) 
(39) 



When CP is conserved, A is a pseudoscalar while 7? 12 are scalars whose mixing is determined 
by the scalar potential. 

Attaching the lepton generation indices, the operator 75, for instance, contains a term 

- -uiu$v 2 cpSpvgv L y, (40) 

which gives neutrino mass after incorporating a coefficient matrix in generations. The phases 
Mi«2 can be removed by redefinition of fields, but will reappear in other terms of T 5 that in- 
volve the Higgs scalars and a lepton pair. These interactions are relatively hard to explore since 
the dominant decays of the scalars generically depend on the details of the underlying theory. 
Furthermore, as we discussed in Introduction, to have any chance at all to discover the mass 
generation mechanism, the mass should be generated from operators of a high enough dimen- 
sion so that the relevant physics scale could be lowered. A promising scenario would be that 
the mass operators are generated at, say, dimension nine, while the lepton-number violating 
operators involving gauge fields are generated at dimension seven by the same physics through 
tree-level or one-loop effects. For these interactions we can say something more certain since 
we know how the gauge bosons interact with the SM fermions. We therefore will concentrate 
on them in what follows. 

Consider first the operators involving gauge field tensors, for instance, M(B). In addition to 
terms involving scalars, it contains the dipole interactions with the Z boson and photon A: 

+ ^UiU* 2 \ 2 cpsp(s w Z^ v -c w A MV )v£<7^ v VLy, (41) 

with cw = cos 0w and sw = sin 6w- A similar dipole term also appears in M 2 (W): 

- -^=uiu* 2 y 2 cpsp(cwZ^ v +swA^ v )v^ x o flv v L y. (42) 

These terms are also contained in the corresponding L operators and barred operators except for 
different factors of wi^, cp and sp. For any given value of /3 there are always operators that are 
not suppressed by its triangular functions. We assign a common coefficient eQ. m . to (the sum) 
of those operators while ignoring factors of order one. While Majorana neutrinos have no dipole 
moments due to CPT invariance, they can accommodate transition dipole moments between 
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different neutrinos IfTTI . Roughly speaking, the upper bounds on the latter are about 10 10 /ls 
or weaker from laboratory experiments lfT2ll and about 10~ 12 /ie from astrophysical arguments 
on energy loss in stars lfT3l . Here /!# = e/ (2m e ) is the Bohr magneton. They translate into a 
bound on the coefficient of the operators: 

10-10 or io~ 12 

C d . m . < 5 ' Le -' Cd - m - ~ ( 6 - 7 OT 31 TeV )" 3 - < 43) 

m e v z 

For the operators involving gauge covariant derivatives, the most interesting interaction is 
the one that contributes to the neutrinoless double beta decay, 

—J 2,4 = / 3 = — mwfaffis^JSxy H , 

K 2A = _ K 3,5 = l„fau 1 ulcpsp£ X y+~-, (44) 

while 7 1 ' 6 , K 1 * 6 , P' 6 do not contain the operator = f^fLyW+^W^ . Here m w = \g 2 \ is the 
W ± boson mass. We assign a common coefficient C xy to (the sum) of these operators. Barring 
exceptional cancellation, their contribution to the subprocess VK~W~ — > ee, ~ C ee m^, should 
not exceed the usual one via the exchange of light active neutrinos, which is experimentally 
bounded and given by ~ m ee /q 2 . Here m ee = Y.j m jU 2 j with mj being the mass of the neutrino 
Vj and U the leptonic mixing matrix, and q ~ (50 ~ 100) MeV is the momentum transfer. The 
upper bounds on m ee lfl4l then imply that 

iC^^L^STeV)- 3 , (45) 
W \ m w 

where we assume for order of magnitude estimation, \m ee \ ~ 0.5 eV and q ~ 100 MeV. 

The operators displayed in the last section also contain other interactions involving multiple 
scalars and gauge bosons, or modify the SM interactions. We leave this more complete phe- 
nomenological analysis for the future work which should better include the effects of multiple- 
fermion operators with a comparable dimension. 
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